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Recent experimental advances in scanning tunneling microscopy make the measurement of the 
conductance spectra of isolated and magnetically coupled atoms on nonmagnetic substrates possible. 
Notably these spectra are characterized by a competition between the Kondo effect and spin-flip 
inelastic electron tunneling. In particular they include Kondo resonances and a logarithmic en- 
hancement of the conductance at voltages corresponding to magnetic excitations, two features that 
cannot be captured by second order perturbation theory in the electron-spin coupling. We have 
now derived a third order analytic expression for the electron-spin self-energy, which can be readily 
used in combination with the non-equilibrium Green's function scheme for electron transport at 
finite bias. We demonstrate that our method is capable of quantitative description the competition 
between Kondo resonances and spin-flip inelastic electron tunneling at a computational cost signif- 
icantly lower than that of other approaches. The examples of Co and Fe on CuN are discussed in 
detail. 



The interaction between conduction electrons and lo- 
calized spins in transition metals with partially filled d 
shells is central to many low-temperature spin effects, 
which may underpin the development of spintronics and 
quantum information devices. When adsorbed on the 
surface of a metallic host, magnetic transition metal 
atoms exhibit various distinctive features in the conduc- 
tance spectrum, which are indicative of many-body scat- 
tering between the conduction electrons and the localised 
spins. These manifest themselves as conductance steps 
at voltages corresponding to the quasi-particle energies 
of specific magnetic excitations and as zero-bias conduc- 
tance peaks, known as Kondo resonances. The first are 
associated to spin-flip inelastic electron tunneling and 
can be described by second order perturbation theory 
in the electron-spin coupling [1], but the second results 
from third order effects due to the electron screening of 
the local spins. 

Recent advances in scanning tunneling microscopy 
(STM) have enabled the detection of many-body scat- 
tering events in Mn [2 , Fe [3 and Co ^ i5j adatoms 
adsorbed on a CuN insulating substrate. The reduced 
symmetry of the surface leads to significant magnetic 
anisotropy, especially for Fe and Co. Fe is found to have 
a large easy-axis anisotropy [D < 0, see Eq. ([2|], lead- 
ing to a ground state spin close to that of the maximum 
z-component of the integer S = 2. This results in four 
evenly spaced conductance steps in the spectrum. For 
Co the large hard-axis anisotropy {D > 0) and the half- 
integer S = 3/2 spin produce a doublet ground state. 
The measured zero-bias Kondo resonance is then due to 
spin transitions between the degenerate ground state lev- 
els. 

Theoretical attempts to reproduce these conductance 
spectra have focused largely on including second order 
scattering events, which cannot account for Kondo res- 
onances but fare well in reproducing the conductance 



steps and their relative intensities [U [6HTTI. Address- 
ing Kondo physics in Co is more involved and one has 
to look for alternative techniques, such as density func- 
tional theory (DFT) informed numerical renormalisation 
group [121 US] • These schemes however are numerically 
expensive. In this letter we extend the perturbative ap- 
proach to the third order in the electron-spin scatter- 
ing and derive an analytic expression for the scattering 
self-energy. This is then implemented within the non- 
equilibrium Green's function (NEGF) formalism [14[ [15] 
for electron transport and used to calculate the STM 
conductance spectra of both Fe and Co on CuN. 

The Hamiltonian describing the STM tip, the mag- 
netic adatom and the nonmagnetic substrate can be di- 
vided into three components, namely a purely electronic 
part, i^e, a purely spin part, i^sp, and an electron- spin 
interaction part, i^e-sp- These are given respectively by 

i^e = ^ Skiaiim + ^0 ^ C^Ca + i^tun , (1) 
kl a 

H,^ = D{Sn' ^ E[{S^f - {Syf] , (2) 

^e-sp = -^sd ^{cl[^]aa^Ca^) ' S . (3) 

a, a' 

In Hq we assume that both the tip (t) and the substrate 
(s) are characterized by a single 5-like band, so that the 
electronic structure (spin-degenerate) is described by the 
creation (annihilation) operators a^i {aki) and the en- 
ergies Ski {I = {t,s}). Similarly the operators cl^/ca 
and the onsite energy eq define the electronic level of 
the adatom. The tunneling between the adatom and the 
tip/substrate is accounted for by i^tun- This term can 
be replaced by an effective energy-dependent self-energy 
i;t,s(£^), which produces a constant broadening of the 
adatom density of states given by rt,s = 7t_a/s-a/^ 
The broadening is thus constant around the Fermi energy 
(E-p) provided that the bandwidth, is large compared 
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to the coupling strength between the two leads and the 
adatom (a), 7t_a/s-a- ^sp contains information on the 
magnetic anisotropy of the adatom and includes both an 
axial (D) and a transverse {E) term, i^e-sp couples the 
conducting electrons to the local spin S through an ex- 
change interaction parameter J^d- 

The calculation of the spin-inelastic conductance spec- 
tra of adatoms adsorbed on surfaces requires the imple- 
mentation of the NEGF scheme including many-body 
scattering effects, which has been rigorously derived in 
reference [1 . In brief, the tunneling current at any ter- 
minal z = t, s and bias voltage V can be calculated from 

/ + 00 
Ti{E) dE , (4) 
-oo 

= ^TTm{E)G<m - p<(^)G>(^)]}, (5) 

where ^>{E) = [1 - /.(^,F)]r, and ^<{E) = 
fi{E,V)Ti. Here fi{E,V) is the bias dependent Fermi 
function in each of the leads, e is the electron charge and 
h the Planck constant. The many body lesser /greater 
Green's functions G^{E) require the calculation of a 
lesser /greater interacting self-energy T^^^{E) which is de- 
rived next. 

Let us consider a nonequilibrium system at finite tem- 
perature described at the level of the Keldysh [16l[ll] for- 



malism. The contour-ordered single particle many-body 
Green's function is defined as 

[G{T,T')U,=-i{\Tc{c,{T)ci,{T')}\). (6) 

This can expanded up to the n-th order \1S^ in the 
electron-spin coupling 

[G(t, t')U, = J2 JdTi...JdT„X 

" c c 

(Q|Tc{i/e-sp(ri) . . . i/e-sp(rn)c.(r)4.(rO}|Q) 
/7(-oo,-oo) ' 

where U is the time-evolution unitary operator and 
the time- averages are now performed over the non- 
interacting ground state |0). 

Exact diagonalization is performed for H^p [Eq. ([2|] 
yielding the spin energies Sm and the states |m) {ftrnn = 
— ^n)- Then by a careful evaluation of the Feyn- 
man integrals in Eq. ^ we arrive at an analytic expres- 
sion for the interacting greater and lesser self-energies, 
^hit^ which includes both contributions to the second and 
the third order in i^e-sp (note that diagrams containing 
fermion loops vanish due to the electron-spin selection 
rules [IE]). This reads (+ corresponds to < and — to >) 



m,n,l i 

- .-(M.) Ee„.(HSin)(n|S^|0«|5V) [H ^ „ ^ I . H ^ J, ^ |] } ■ (8) 

I 



where /cb is the Boltzmann constant, T is the temper- 
ature and Pi is the population of the |/) spin state. If 
we now assume that the adatom is much more strongly 
coupled to the substrate than to the STM tip (75 ^ 7t), 
we can approximate its density of states around Ep with 
a constant, p = (rs/27r)/(6:o + r^/4), where Sq > E^. 
The weak coupling to the STM tip also ensures that 
the spin system remains always close to equilibrium \10\ . 
This means that the adatom state always resides close to 
the ground state, i.e. that Po^l- The matrix elements 
{m\S'^\n), with i = {x^y^z}, determine the intensity of 
a given transition between an initial state n and a final 
state m. We also note that the strength of the 3rd order 
interaction depends on the dimensionless parameter pj^d 
in agreement with the expression of Ref. [19j. 

Having found a close formula for the interacting self- 
energy we can now proceed with the full NEGF formal- 
ism [1 and calculate the conductance spectra of various 



adatoms on insulating substrates. In order to explore the 
Kondo physics contained in the 3rd order contributions 
to S^^, we consider both Co and Fe adsorbed on CuN, 
for which we can compare with experiments [3]-[5]. 

Figures [1] and [2] show the calculated conductance spec- 
tra respectively for single Co and Fe atoms on CuN as 
the dimensionless parameter a = pJsd is varied (note that 
the spectra are normalized relatively to the elastic only 
conductance. Go)- Many of the parameters of the model 
used in the present simulations have been extracted from 
the DFT calculations [T2|- The broadening due to the 
coupling to the substrate (we neglect that due to STM 
tip) is set at Fg = 0.1 eV for both Co and Fe. In order to 
ensure consistency in our approximations (Fg = 7s_a/^) 
we also set 7s-a = 1.5 eV and the substrate bandwidth 
to = 20 eV. We then choose eo = I eV to fulfin the 
criterion sq ^ E^ = 0. The magnitude of Jgd for both Fe 
and Co is held constant at 0.5 eV (note that Jgd =0.5 eV 
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FIG. 1: Normalised conductance spectrum for Co (S = 3/2) on 
CuN as q; = pJsd is increased. Note the emergence of a Kondo 
resonance at zero-bias, i.e. at the Fermi level. 



FIG. 2: Normalised conductance spectrum for Fe (S = 2) on CuN 
as q; = pJsd is increased. Note the transition from a constant 
conductance at the inelastic steps to a conductance logarithmic 
decaying. 



is comparable to what calculated from DFT for a Ni im- 
purity sandwiched between Au leads [20]). Finally the 
empirical parameters describing the axial, and trans- 
lational, E anisotropics are taken from Refs. [3 ^ and 
are I^co = 2.75 meV, Eqo = meV, D^e = —1.53 meV 
and Epe = 0.31 meV, while the adsorbed atoms spins are 
5'co 3/2 and S'pe = 2. 

The parameter a is varied by altering the coupling of 
the adatom to the substrate, 7s- a, thus by changing the 
value of the density of states at the Fermi level, p. A full 
diagonalisation of H^p for Co gives a set of four (2S'co + l) 
eigenvalues and eigenvectors. In particular the presence 
of hard-axis anisotropy results in the following energy 
manifold = {0.69,0.69,6.19,6.19} meV, i.e. in a 
doubly degenerate ground state. It is then found that 
transitions between the degenerate ground states only 
become allowed on inclusion of the third order term in 
Eq. ([8| because of the selection rules imposed by the the- 
ory through the matrix elements {m\S'^\n). Such a tran- 
sition appears in the spectrum of figure [l] in the form of 
a zero-bias Kondo peak as the value of a is increased. 
The greater is the density of states of the adatom, the 
greater becomes the third order Kondo contribution to 
the conductance. This also produces a logarithmic rise in 
conductance both for the ground and for the first excited 
state transitions (detected at around 6 mV). These fea- 
tures matche closely what is observed in the experiments 
of Refs. |4| |5j. A notable difference is that our calcu- 
lation seem to show a peak considerably smaller than 
that in experiments. Such a discrepancy may arise from 
two-body electronic (Hubbard-like) contributions to the 
Hamiltonian, which we neglect here. 

Results for Fe are presented in Fig. [2] This 
time the five eigenvalues of H^p are = 
{-6.30,-6.12,-2.46,-0.60,0.18} meV, so that the 



ground state is non-degenerate. At zero magnetic field 
all transitions available by the 3rd order expansion are 
resolved. As a is increased the 3rd order Kondo term 
produces a logarithmic peak at each of the allowed 
conductance steps. This feature agrees well with exper- 
iments and it was previously explained by invoking a 
non-equilibrium population of the Fe spin states [10]. 
This explanation however confiicts with the fact that the 
same effect is seen also for low currents [3^ , namely when 
the adatom spin state resides close to the equilibrium. 
Therefore, we re-interpret the experimental data as 
the manifestation of the 3rd order Kondo effects in the 
spin-fiip inelastic tunneling spectra. 

As a final test for our 3rd order self-energy we consider 
a Co-Fe dimer with exchange coupling between the ions, a 
situation already investigated experimentally in Ref. [5]. 
In our formalism this translates in including an additional 
Heisenberg-like term to H^p 

= Jdd(Sco • Spe), (9) 

where Jdd is the exchange coupling constant. Figures [3] 
and [4] show the calculated conductance spectra when the 
tip is positioned respectively over Co or Fe. Note that in 
general Jdd is small so that any changes in the electronic 
levels of the combined Co-Fe system are not resolved in 
the mV range and can be neglected [l]). In general we no- 
tice a change in the conductance spectrum of each atom 
as they are brought closer together, i.e. as Jdd increases. 
The conductance step at 0.18 mV for Fe decreases with 
increasing Jdd, while Fe itself acts as an effective mag- 
netic field that splits the zero-bias Kondo resonance of 
the Co spectrum. These effects are both observed in the 
experiments of Ref. [5 . Notably Fe does not simply act 
as a source of magnetic field on the Co, as seen in the in- 
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FIG. 3: Normalised conductance spectrum of Co when it is ex- 
change coupled to Fe on CuN as J^a is increased. The Co Kondo 
peak splits as the Fe atom acts as an effective magnetic field. The 
insert zooms in the zero-bias region. 

set of figure [3] for J^d = 0.03 meV. In fact one can clearly 
observe an additional Kondo peak emerging at zero-bias 
in between the principally split peaks. This is a unique 
feature of the exchange coupling between Co and Fe. In 
fact there are now (2S'co + l) x (26'Fe + l) = 20 eigenvalues 
and additional allowed transitions appear at each atomic 
site. For instance for large Jdd the zero-bias region be- 
comes completely dominated by a conductance dip. This 
originates from the opening of a spin transition between 
the ground state at -5.686 meV and the first excited state 
at -5.379 meV. Such a transition has a spectral intensity 
much larger than that of the Kondo resonance, which 
therefore disappears from the spectrum. 

In conclusion we have studied the effects of including 
3rd order contributions in the electron-spin coupling to 
the conductance spectra of transition metal atoms ad- 
sorbed on a CuN substrate. In particular we have derived 
a close expression for the 3rd order electron-spin self- 
energy within the NEGF formalism and a single band 
tight-binding model incorporating local Heisenberg ex- 
change to quantum spins. Two main features in the 
conductance spectra arise from our formalism, namely 
a logarithmic decay of the conductance as a function of 
bias and zero-bias Kondo resonances. We obtain a good 
agreement with available experimental data for both iso- 
lated Co and Fe atoms and for an exchanged coupled 
Co-Fe dimer. Importantly the low computational effort 
needed by our method makes it a valuable alternative 
to full many-body treatments in describing spin inelastic 
phenomena at the atomic level. We also believe that our 
proposed scheme is amenable to be combined with first 
principles methods, i.e. it can form the basis for a fully 
quantitative theory of spin scattering in nanostructures. 

This work is sponsored by the Irish Research Council 
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FIG. 4: Normalised conductance spectrum of Fe when it is ex- 
change coupled to Co on CuN as J^d is increased. Note that the 
intensity of the conductance step at 0.18 mV decreases with in- 
creased exchange coupling. 
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